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ABSTRACT 

5? 



Context. In the single degenerate (SD) scenario of type la supernovae (SNe la) the non-degenerate companion star survives the 
supernova (SN) explosion and thus should be visible near the center of the SN remnant and may show some unusual features. 
Therefore, a promising approach to test progenitor models of SNe la is to search for the companion star in historical SN remnants. 
Aims. Here we present the results of three-dimensional (3D) hydrodynamics simulations of the interaction between the SN la blast 
O { \ wave and a MS companion taking into consideration its orbital motion and spin. The primary goal of this work is to investigate the 

I , rotation of surviving companion stars after SN la explosions in the WD+MS scenario. 

Methods. We use Eggleton's stellar evolution code including the optically thick accretion wind model to obtain realistic models of 
companion stars. The impact of the supernova blast wave on these companion stars is followed in 3D hydrodynamic simulations 
| employing the smoothed particle hydrodynamics (SPH) code GADGET3. 

Results. We find that the rotation of the companion star does not significantly affect the amount of stripped mass and the kick velocity 
caused by the SN impact. However, in our simulations, the rotational velocity of the companion star is significantly reduced to about 
14% to 32% of its pre-explosion value due to the expansion of the companion and the fact that 55% - 89% of the initial angular 
momentum is carried away by the stripped matter. 

Conclusions. Compared with the observed rotational velocity of the presumed companion star of Tycho's supernova, Tycho G, of 
■ ~ 6 km s _1 the final rotational velocity we obtain in our simulations is still higher by at least a factor of two. Whether or no this 

\^ ■ difference is significant, and may cast doubts on the suggestion that Tycho G is the companion of SN 1572, has to be investigated 

| in future studies. Based on binary population synthesis results we present, for the first time, the expected distribution of rotational 

velocities of companion stars after the SN explosion which may provide useful information for the identification of the surviving 
companion in observational searches in other historical SN remnants. 

. Key words, stars: rotation- supernovae: general-hydrodynamics-binaries: close 

- 1— . 1 . Introduction In principle, there are various possibilitie s for the evolution 
^ . towards an explosion (see lWang & Hanll2012l for a review). One 
J-H ■ option is a merger of two CO WDs with a combined mass in ex- 
Type la SNe (SNe la) are used as cosmic distance indicators cess of the Chandrasekhar mass, which may explode as a SN la , 
since their luminosity can be calibrated based on the empir- the 'double de generate' (DD) scenario (Ibe n&Tutukovi ri984: 
ical relation between light curve shape and peak luminosity IWebbinklll984l) . The DD model can explain in a natural way the 
dPhillipsI 1 1 993l: iPhillips et al.lll999h . This has provided the first lack of hydrogen in SNe la. However, only a few DD systems 
evidence for the a c celerating expansion of the present univers e have been found whose orbital period is short enough that they 
dRiess et al.l 1 19981: iPerlmutter et al.l 119991: iLeibundgutl 120081) . will merge within a Hubble time, but in none of them the com- 
SNe la are also believed to be important contributors to the cos- bined mass exceeds the Chandrasekhar ma ss limit. On the other 
mic nucleosynthesis and they are sources of kinetic energy in hand, in recent numerical simulations of iPakmor et al.l (1201 Ol 
evolution process of galaxies. Although a detailed understand- 1201 ll) it was found that the violent merger of a pair of white 
ing of the nature of thei r progenitors and the physics o f explo- dwarfs with equal masses of ~ 0.9 M© can directly trigger a 
sion is still lacking (see iHillebrandt & Niemeyei]|2000l for a re- thermonuclear explosion that resembles sub-luminous 1991bg- 
view), there is consensus that SNe la arise from thermonuclear like SNe la. Moreover, it was shown that the violent merger 
explosio ns of carbon/oxygen white dwarfs (CO WD s) in binary of two CO WDs with masses of 0.9 M© and 1.1 M© produces 
systems dHoyle & Fowledll960tlNomoto et al.lll997l) . lightcurves and spectra which are in good agreement with those 
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of normal SNe la (iPakmor et al.ll2012bl) . supporting the DD sce- 
nario. 

On the other hand, a rather massive WD may accrete 
hydrogen-rich material from a non-degenerate binary com- 
panion until it approaches the Chandrasekhar mass, the 'sin- 
gle degenerate' (SD) scenario. In this case, the binary com- 
panion could be a main- sequence (MS) star (WD+MS chan- 
nel), a slightly evolved sub giant star or a red-gia nt (RG) star 
(lHan & Podsiadlowskil[20Q4tlWang & Hanll2010bl) . The lack of 
hydrogen in observed SN la spectra can be seen as trouble- 
some for the SD scenario since the companions are hydrogen- 
rich stars. With more realisti c MS companion star models than 
those used in previous work, iLiu et al.l (1201 2l) performed three- 
dimensional (3D) smooth particle hydrodynamics (SPH) sim- 
ulations of the interaction between SN la ejecta and the MS 
companion star. They found that in all cases they considered 
more than 0.1M© was stripped from the hydrogen-rich compan- 
ion and was mixed into the SN ejecta, which is in disagree- 
ment with most recent observational constraints on the pres- 
ence of hydrogen in SN la from nebular spectra ~ 0.01 M 
(lLeonardl2007l:IShappee et al1l2013l) . Moreover, no other similar 
hydrodynamics simulations showed an am ount of stripped mass 
below this strong observational limit (see iMarietta et al.l [2 000 ; 



Pakmor et al. 2008; Pan et al. 20101 iRicker et alJ2010l:lPan et al] 



2012bl) . On the other hand, there is evidence from observations 
that the progenitors of at least some SNe la come from the SD 
channel. For instance, evidence for the presence of circumstellar 
matter and features indicative for an interactio n between the SN 
and circumstellar matter were found recently dPatat et al.l 20071; 



ISternberg et aHl201lHFolev et al.ll2012HDildav et alfco'm 

The merger of two WDs leaves no remnant after a SN la 
explosion. In contrast, in the SD scenario the companion star 
survives the explosion and, in principle, can be identified due to 
its peculiar spatial velocity, its rotation , effective temperature , 
luminosity or composition (iHanl 120081: IWang & Hanfl2010al) . 
Therefore, it is a promising approach to test the progenitors of 
SNe la by directly searching for the surviving companion star in 
galactic SN remnants (SNRs). There are several ways by which 
the SN blast wave modifies the properties of the companion. 
First, the SN strips off matter from the surface of the compan- 
ion and injects thermal energy into it during the interaction. This 
causes the surviving companion to expand and lowers its surface 
gravity. Secondly, after the impact the companion star's surface 
will be enriched with heavy elements (e.g., Ni, Fe or Ca) from 
the inner part of the SN ejecta which should show up in its sp ec- 
trum (see iGonzalez Hernandez et aL 

Finally, after the explosion the companion star retains its pre- 
explosion orbital velocity, giving it a peculiar velocity compared 
to other stars in the vicinity. 

Tycho Brahe's SN 157 2 isaS N la that exploded in the Milky 
Way. iRuiz-Lapuente et al.l (120041) have analyzed the stars within 
a circle of 0.65 arcmin radius of the center of the SNR up to 
an apparent visual magnitude V = 22. They found a star, Tycho 
G, similar to the Sun in surface temperature and luminosity but 
with a lower surface gravity than a MS star. It has a significant 
peculiar velocity in radial and proper motion, and moves at more 
than three times the mean velocity of the other stars in the field. 
Therefore, they suggested that Tycho G star could be the sur- 
viving companion star of SN 1572. However, since then it has 
bee n noted that Tych o G does not show any spectral peculiari- 
ties (llhara e t al. 2007 1) and t hat it is apparently not out of thermal 
equilibrium (lHowellll201ll) . iFuhrmannl (120051) also claimed that 
it might be a thick-disk star coinciden tally passing in the vicinity 
of the remnant of SN 1572. Recently jGonzalez Hernandez et al.l 



(l2009h found that Tycho G has an overabundance of Ni rela- 
tive to normal metal-rich stars, and they suggest that Tycho G 
could have captured the low- velocity tail of the SN 1572 ejecta, 
which upholds Tycho G as a surviving companion star again. 
H owever, the measured [N i/Fel ratio from a more recent study 
of iKerzendorf et al] (120121) seems to be not so unusual with re- 
spect to field stars with the same metallicity. 

In the SD case, because of the strong tidal coupling of a 
Roche-lobe filling donor, the donor star rotation is expected to be 
tidally locked to its orbital motion. This forces the binary donor 
star to have a sp in corresponding to the orbita l frequency of the 
binary system (lHan & Podsiadlowskil 120041; IKerzendorf et al] 
120091) . After the SN explosion, the companion is released from 
its orbit and continues to rotate. Therefore, fast post-explosion 
rotation might be a signature of the donor star. The survivor 
coul d be in rapid rotation which would be measurable easily 
(see IKerzendorf et al.ll2009l l2012h . With the HIRES instrument 
on the Keck-I telescope and with Sub aru high-resolutio n spec- 
trosco py of star G in the Tycho SNR IKerzendorf et al.l (120091 
120121) measured the rotational velocity of Tycho G and ob- 
tained only -6+1.5 km s" 1 . Thus, they concluded that Tycho 
G is unlikely to be the surviving companion star of SN 1572 
because it does not have the expected high rotational velocity 
(IKerzendorf et al.ll2009h lHowelll[20TTh . However, the impact of 
the SN ejecta strips off material from the surface of the compan- 
ion star, which can significantly reduce its angular momentum, 
thereby lowering the rotational velocity. Moreover, the SN im- 
pact and heating bloat th e companion star, which will also slow 
down its spin. Recently, iPan et al] (1201 2al) computed the post- 
impact evolution of the remnant star produced in their multi- 
dimensional adaptive mesh refine ment simulations of WD+MS- 
like models with the FLASH code (IPan et al.l2012bl) . In line with 
the arguments given above, they conclude that Tycho G cannot 
be eliminated as being a promising proge nitor candidate, b ased 
on its low rotational surface velocity only (P an et al ] l2012al) . 

In this paper, we present the results of 3D SPH simulations 
of the impact of SN la ejecta on a MS companion star includ- 
ing the orbital motion and spin of the companion. The purpose 
of the work is to investigate how the SN la impact changes the 
rotation rate of the companion which will be useful for identi- 
fying the surviving companion star in SNRs. The paper is orga- 
nized as follows. In Section [2l we describe the codes used and 
show an example of the initial companion star models used in 
our simulations. In Section [3] the effect of rotation on the in- 
teraction between SN la ejecta and its MS companion are dis- 
cussed. Post-impact rotation features of a MS star are shown 
in detail in Section |4] Furthermore, the results of our simula- 
tions are compared with the observed rotation of star Tycho G 
in Section [5] Next in Section O we discuss the distribution of 
the rotational velocities of the remnant stars after thermal equi- 
librium is reestablished. Finally, in Section [7] we summarize our 
results. 



2. Numerical method and model 

2.1. Numerical codes and initial setup 

In order to obtain consistent MS comp anion star mode ls at 
the onset of the SN explosion, similar to ILiu et al.l (I2012I) , we 
used Eggleton's one-dimens ional (ID) stellar evolution code 
(lEggletonlll97li I1972L 1 1973b to follow the binary evolution of 
a SD progenitor system in detail. Roche- lobe overflow (RLOF) 
is trea ted in the code as described by lHan & Podsiadlowskil 
( 2004). The opacity tables used in our calculations are com- 
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piled bvlChen & Toutl (f2007l) f rom llgiesias & Rogers! (1 19961) and 
I Alexander & Ferguson (1 1994b . We use a typical Population I 



composition with hydrogen abundance X = 0.70, helium abun- 
dance Y = 0.28, and metallicity Z = 0.02. We set a = l/H F , the 
ratio of mixing length to the local pressure scale heig ht, to 2 and 
set the convective overshoot ing parameter S 0Y to 0.12 dPols et alJ 
119971: ISchroder et al.lll997h . which roughly corresponds to an 
overshooting length of ~ 0.25 pressure scale height. 

Furthermore, the lat e st version of the SPH code GADGET- 
3 (ISpringel et al.l 120011: ISpringell 120051) is employed to sim- 
ulate the impact of SN la explosions on their MS compan- 
ions by including the orbital and spin velocities of companion 
stars. The GADGET code was originally used for cosmologi- 
cal simulations, but it has been modi fied to make it appli cable 
to stellar astrophysics problems (see iPakmor et al.l 1201 2al) and 
it has been used successfully to capture the main dynamical ef- 
fects of the SN impac t on its companion star (see lPakmor et al.l 
20081: iLiu et aHl2012h . In our simulations, we aim at determin- 
ing the rotational velocity of the companion star after the impact. 
Therefore the fact that no matter leaves the computational do- 
main and that momentum, energy and angular momentum are 
strictly conserved are a crucial advantage of SPH over grid- 
based methods for the problem under investigation. 

In this work, th e basic setup for the GADGET code is al- 
most the same as in ILiu et alJ (120121) . The smoothing length is 
chosen such that a sphere of its radius encloses 60 neighbor- 
ing particles. All particles have the same mass. The gravitational 
softening length is equal to the smoothing length. To reduce 
numerical noise introduced by mapping the initial model, the 
MS companion star is relaxed for 1.0 x 10 4 s (several dynami- 
cal timescales) be fore we start the actu al impact simulation. The 
W7 SN la model dNomoto et al.lll984l) is used for the SN explo- 
sion, and the orbital separation is taken from the ID consistent 
binary-evolution calculations. 

Additionally, in this work, we assume that the companion 
star co-rotates with its orbit due to strong tidal interaction. Thus 
the spin period of the companion star is locked with its orbital 
period. Furthermore, the orbital and spin velocities of the com- 
panion stars are included in our impact simulations in order to 
study the post-impact rotation rate of surviving star. 

The system of coordinates is chosen as follows. We set the 
x - y plane to be the orbital plane of the binary system and as- 
sume a circular orbit. The z-axis is chosen as the spin axis of 
the companion, and the positive z-axis is the direction of the 
angular momentum. Finally, we assume that the companion star 
rotates as a rigid body at the moment of the SN explosion. 

2.2. Progenitor model 



As in ILiu et al.l (12012b (see also lHan & Podsiadlowskil 120041: 
IWang & Hanll2010al) . we start our ID consistent binary calcula- 
tion when the WD+MS system has been formed, the mass trans- 
fer then occurs through RLOF when the companion star fills 
its Roche lobe. Instead of solving the stellar equations for the 
WD star when we trace the detailed evo lution of the compan ion 
star, the optically thick wind model of lHachisu et al.l (1 19991) is 
adopted. 

Figure [T] shows the evolution of a binary system with an 
initial donor-star mass of M l 2 = 2.0 M , an initial mass of the 
CO WD of M^ D = 0.8 M© and an initial orbital period of 
log(PVday) = 0.2. In such a system, the companion star has 
an orbital velocity of 138 km s" 1 at the moment of the SN la ex- 
plosion. Assuming that the rotation of companion is locked with 
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Fig. 1. Binary-evolution calculation of the MS_110 progenitor 
model. In panel(a), the solid, dashed and dash-dotted curves 
show the mass-transfer rate, M2, the mass-growth rate of the 
CO WD, M co , the mass of the CO WD, M WD , respectively. In 
panel (b), the evolutionary track of the donor star is shown as a 
solid curve and the evolution of the orbital period is shown as 
a dashed-dotted curve. The dotted vertical lines in both panels 
and asterisks in panel (b) indicate the position where the WD 
explodes as a SN la. 



the orbital motion due to the strong tidal interaction, the rota- 
tional velocity, v rot , and its orbital velocity, Vprb.2, o bey a simple 
relation as follows (see also Kerzendorf et al. 2009): 



Mi + M 2 . , , 

Vrot = H JWVorb,2, 

Ml 



(1) 



where q = M2 /M\ is the mass ratio of the binary system at the 
moment of the explosion, and f(q) is the ratio of t he Roche- 
lobe r adius of companion star to the orbital separation (Eggleton 
Il973l) . Therefore, we calculate the rotational velocity of the 
companion star to be ~ 110 km s" 1 at the time of the ex- 
plosion. 

In this model, because hydrogen burning on the WD is un- 
stable before the SN explosi on, the system may be observed as 
a U S co-type recurrent nova (lHachisu et a l. 2008; Meng & Yangl 
12010b . With our consistent binary-evolution calculations, we ob- 
tain a CO WD+MS progenitor that is similar to the U Sco binary 
system with a mass of the secondary of 1 . 1 8 Mq and an orbital 
period ~ 1 day (see Table [T]). This model is named "MS_110", 
the "MS" indicating a CO WD+MS binary system, the "110" 
means the rotational velocity of the companion star of ~ 
110 km s" 1 at the moment of the explosion. U Scorpii is one of 
the best-observed recurrent novae, and it has been suggested as a 
progenitor of a SN la because its white dwarf mass is close to the 
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Fig. 2. Density distribution of the SN ejecta and the companion 
shows the logarithm of the density. The plots are made using the 



Chandrasekhar mass (lHachisu et al.ll200(A iThoroughgood et all 
120011; |Podsiadlowskill2003l) . In order to investigate th e physical 
parameters of the recurrent nova, lHachisu et al.1(l2000l) have suc- 
cessfully modeled the theoretical light curve for the outburst of 
the U Scorpii system with a model consisting of a WD mass of 
Mwd = 1.37M©, a secondary mass of M2 = 0.8 - 2.0M© and an 
orbital inclination of ~ 80°. 

3. Effects of rotation 

Based on the system similar to U Sco (MS_110) that was dis- 
cussed in Section 12.21 we performed a 3D SPH impact simula- 
tion which included the spin velocity of = HOkms" 1 , and 
the orbital velocity of v^ 2 = 138 km s" 1 to study the effect of 
asymmetry due to the orbital motion and spin of the companion 
star. We used a total of 2 x 10 7 particles to represent the compan- 



rotation of surviving companions 




star of a non-rotating version of model M S_110. The color scale 
freely available SPLASH tool dPricdl2007b . 



ion star only (which corresponds to a total number of particles 
in the simulation of ~ 4 x 10 7 ). All particles had the same mass 
of 5.9 x 10" 8 Mp dot- The SN propertie s were taken from the W7 
explosion model dNomoto et al .11 19841) with an initial orbital sep- 
aration of 3.75 x 10 11 cm. 

The temporal evolution of the spatial density distribution of, 
both, the companion and the SN ejecta MS_1 10 model with and 
without including the orbital motion and the rotational velocity 
of the companion are shown in Fig. [2] and Fig. [3] respectively. 
Some small morphological differences between the rotating and 
non-rotating model are seen due to the symmetry-breaking ef- 
fects of orbital motion and rotation of the companions star (see 
Fig. E]). 

At the end of the simulation, a mass of ~ 0.23 M© (~ 19% of 
the total companion mass) is stripped from the companion star 
in the MS_110 model (see Fig. [4]). The orbital motion and spin 
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Fig. 3. Same as Fig.[2l but including the orbital motion and the spin of the companion star. The color scale shows the logarithm o f 
the density. Some differences become visible at late times. The plots are made using the freely available S PLAS H tool ( Price l2007l) . 



of the companion star only leads to a ~ 2% larger unbound mass 
(see Fig.©, but basically the same kick velo city (~ 60kms~ 1 ) . 
The result differs from what was reported by iPan et all (120 1 2bl) . 
who found that 16% more mass can be stripped if the orbital 
motion and spin is included. In Section [5j again, we calculate 
the amount of unbound mass for the model MS _ 160 which has a 
higher rotational velocity of 160 km s" 1 , but still only a 4% dif- 
ference is found. Moreover, we do not find that the rotation of 
the companion star significantly affects the post-impact velocity 
distribution of unbound material and SN la ejecta (see Fig. [5}. 
Most of the stripped hydrogen-rich material is confined and hid- 
den close to the center of the SN ejecta. The detection of a hydro- 
gen line may be possible only when the high velocity SN ejecta 
becomes transparent. 

To summarize, the main results (e.g., the unbo und mass, the 
kick velocity) are similar to our previous work (see Pa kmor et al.l 



120081: iLiu etai1l2012h even if the orbital and spin velocity of the 
secondary are taken into account. This is not surprising, since 
the orbital and spin velocity of the companion star are obviously 
much lower than the typical values of the expansion velocity of 
the SN la ejecta (~ 10 4 km s" 1 ). Therefore, rotation cannot affect 
the basic physics of interactions between the SN ejecta and the 
companion star significantly. 



4. Post-impact rotation 

In this section, we investigate how the SN la explosion affects the 
rotation of the companion star in our SPH impact simulation s by 
using again the MS_1 10 model discussed in detail in Section 
as a typical case. 
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4.1. Temporal evolution of the overall rotation 

At the time of the SN explosion, the companion star is expected 
to rotate rapidly at the same frequency as the binary system due 
to the strong tidal interaction during the RLOF phase. After the 
SN impact, however, it is uncertain whether it is still possible to 
use signatures of rapid rotation to single out candidates for donor 
stars. This depends on the total amount of angular-momentum 
taken by the stripped material and extreme expansion of the 
donor star due to SN ejecta heating. 

By including the orbital motion and spin of the MS_110 
model, we investigate the total angular momentum that can be 
carried away by the stripped mass, studying the post-impact evo- 
lution of rotational velocity of the surviving companion star. 
During the interaction, the collision of the SN la ejecta brings 
the companion star out of thermal equilibrium. Simultaneously, 
the spherical symmetry of the companion star is also strongly 
affected, deforming the shape of the star (the details of density 
evolution of a remnant star are shown in Section |4~2] ». Therefore, 
it is difficult to determine the real surface of a surviving com- 
panion star and to obtain its overall rotational velocity after the 
SN impact. For this purpose, the two following steps are done 
to estimate the surface of the surviving companion star in this 
work. 

a) We divide the surviving star into several hundred spherical 
shells. The density of the SPH particles in each shell are aver- 
aged to calculate a value for that shell. Next, if the density fluc- 
tuation in a shell is too large, it is ignored because the very outer 
shells are very poorly resolved anyway (SPH particles in these 
shells are too rare to reasonably resolve the structure of the star). 
The density profiles of the MS_1 10 model as a function of explo- 
sion time are shown in Fig. [6] At this stage, we obtain the surface 
of surviving star at R\, where the tangential velocity with regard 
to the initial rotation axis, z-axis, is v r ° l . 

b) In Fig. [61 a sharp decrease in density is seen at the outer 
layer of the surviving companion star. Here, we simply choose 
the position R2 (~ 95% of R\) of the sharp density jump as the 
real surface of the star (see the vertical dotted line in Fig. O. 
Furthermore, we take the rotation velocity, v™ 1 , at this surface, 
R2, to denote the overall apparent rotational velocity of the sur- 
viving companion star, v rot (i.e., equatorial rotational velocity). 
In the following Section 14.21 it is shown that the rotational ve- 
locity is not changing significantly in the outer parts of the stel- 
lar envelope. The difference between v r ^ and v r ^ is < 5%. On 
the other hand, it is also difficult to determine the exact position 
of the photosphere of the surviving companion star. However, 
Figure [6] shows that the density at the position R2 is always 
above the lower density-limit of the photosphere of the sun 
(~ 10" 7 gem -3 ). Therefore, using the rotational velocity of the 
star, v™ 1 , to estimate its observed rotational velocity should be a 
good approximation. 

Figure |7] shows the temporal evolution of / spm , /, co, v rot and 
Abound of the companion star in the impact simulation for the 
MS_110 model. Here, J spm , /, oj, v rot and M boim d are the total 
angular momentum, the total moment of inertia, the angular ve- 
locity and the rotational velocity at the equator, and the bound 
mass of the companion star, respectively. In the MS_110 model, 
the companion star has an initial total angular momentum of 
2.25 x 10 50 gcm 2 s -1 at the moment of the explosion. At the 
end of the simulation, 0.23 M (~ 19% of the initial mass) are 
stripped-off due to the SN impact, carrying a spin angular mo- 
mentum of 1.63 x 10 50 g cm 2 s _1 (~ 72% of that before the SN 
explosion). Therefore, the surface rotational velocity of the star 



0.25 




Fig. 4. Unbound mass from time of explosion with (solid line) 
and without (dotted line)the orbital rotation and spin of of the 
companion star for model MS_1 10. 



is significantly reduced to ~ 25 km s 1 from an initial rotational 
velocity of ~ 1 10 km s" 1 . 

The SN impact takes away both mass and angular momen- 
tum, which causes the total angular momentum of the compan- 
ion star to decrease sharply. Therefore also the angular veloc- 
ity and the rotation velocity (see the dashed lines in Fig. [7]) de- 
crease significantly. Furthermore, the total angular momentum 
stops decreasing (see dashed line) as the amount of the stripped 
mass (solid line) reaches a constant value after 5000 s, which 
was shown in Fig.[7k. As the companion star is puffed up due to 
energy deposition from the impact of the SN ejecta, its outer lay- 
ers expand, remarkably increasing the moment of inertia of the 
star (see Fig. (TJ^). This explains why the angular velocity keeps 
slowly decreasing after 5000 s (see dashed line in Fig. |7t) al- 
though there is no additional mass-loss and angular-momentum- 
loss at this moment. However, this expansion does not reduce 
the rotational velocity of companion significantly as is shown by 
dashed line in Fig.|7Ji. 



4.2. Post-impact radial-velocity distribution in the companion 
star 

In our simulations, the initial rotation of the companion star was 
set up as a rigid-body rotation, the rotation axis being the z-axis. 
At the moment of the SN explosion, the companion star is spher- 
ically symmetric, having the same angular velocity from the cen- 
ter to the surface. Therefore, the rotational velocity increases lin- 
early with radius (see right diagram on the top row of Fig. |8}. 

Figure [8] shows the density distribution (left column) and the 
radial profiles of the angular velocity (solid curve), oj, and the ro- 
tational velocity (dashed curve), v rot , of the surviving companion 
star at different times after the explosion for the MS_1 10 model. 
After the SN impact, heating by the SN ejecta puffs up the com- 
panion star, causing its envelope to expand considerably. This 
produces a object with a compact inner core and a low-density 
outer layer (see left column in Fig.[8j. The star starts to relax to- 
wards a spherical state about 2 x 10 4 s after the impact. However, 
the density in the outer layers varies by as much as a factor of 
~ 3 in different directions (see bottom r ow in Fig. [8), which 
was also seen in previous simulatio ns (see iMarietta et al.l l2000: 
iPakmor et aDl2QQ8fc iLiu et al.ll2012h . At this moment, however, 
the sound crossing time is still longer than the time we follow in 
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Fig. 5. Velocity distribution of stripped material that originally belonged to the companion star (left figure) and SN ejecta (right 
figure) for the non-rotating (solid lines) and rotating (dashed lines) MS_110 model. 



the simulations (~ 5.6 hours). Therefore, the star does not reach 
a spherical state. 

The right column of Fig. [8] shows that the equatorial rota- 
tional velocity of the companion star drops to ~ 25 km s" 1 from 
the original rotational velocity of 1 10 km s" 1 due to the extreme 
expansion of the star and the angular momentum losses. The an- 
gular velocity near the center of the star also decreases below its 
initial value of 0.8 x 10" 4 rads _1 due to the effect of the shock 
running through the center, but it becomes stable at late times 
(see Fig. [9]). Additionally, Fig. [8] shows how the radial rotation 
profiles of the companion star are affected during the collision 
with the SN la ejecta (see also Fig. |9j). After the SN impact, the 
companion star is no longer in rigid-body rotation and its outer 
layers exhibit some features of differential rotation. At the equa- 
tor the companion star rotates at a different angular velocity than 
at higher latitudes. The detailed radial distribution of various ro- 
tational properties of model MS_110 2 x 10 4 s after the explo- 
sion are shown in Fig. [TUl The angular velocity decreases with 
increasing latitude of the surviving companion star: the equator 
rotates with a period of ~ 14 days, near the poles it is as small 
as 2 days. In order ensure that changes of the rotational profile 
of the companion star are purely caused by the impact of the SN 
explosion, we have also done a test run without the SN. It shows 
that without the SN the rotational velocity profile of a compan- 
ion star does not change at all over the time of the simulation 
(2 x 10 4 s). In our simulation, we use the rotational velocity at 
the surface of the post-impact remnant star as the potentially ob- 
servable velocity. If the photosphere of star should be located at 
a smaller radius rather than the surface, the rotational velocity 
would be a little bit higher than ~ 25 kms -1 . However, the dif- 
ferences of the rotational velocity are not significant for a large 
range of outer layers (see Fig. [8]and Fig. [9]>. Also, in this work, 
the initial rotation of the companion star is set up in rigid rota- 
tion, which might be inconsistent with the realistic rotation of a 
star. A more realistic rotation profile of the star is required to re- 
place the rigid rotation used here in forthcoming investigations, 
although we do not expect this to change our main results. 




r=R 1 
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© © 



^© " 



2.0 3.0 4.0 5.0 

Radius [ 10 11 cm 1 
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Fig. 6. Density profile of the MS_110 model for different times 
after the SN explosion. The black solid line corresponds to the 
initial density profile of the star at the moment of the explosion. 
All black dashed lines from (a) to © show the radial density dis- 
tribution of the star (0.4, 0.8 •• • 2.0) x 10 4 s after the impact in 
time intervals of 4000 s. 



Figure [TT] shows the temporal evolution of the radial distri- 
bution of the angular momentum of the companion star in model 
MS_110. The different colors belong to different times since the 
SN explosion. Each colored solid line shows the details of the 
radial profile of the angular momentum at that given time. The 
angular momentum, J r , was calculated by summing up the total 
angular momentum of all bound particles interior to the corre- 
sponding cylindrical surface at radius r. The dashed-dotted line 
shows the change in total angular momentum as time progresses. 
Most of the angular momentum of the surviving companion star 
is concentrated in a small inner core. Outside this core, how- 
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Fig. 7. Temporal evolution of some properties of the companion star after the SN impact in MS_110 model. The black solid line 
corresponds to the total bound mass of the companion star, Mb oim d- The dashed curves show the temporal changes of the total 
angular momentum (J spm ), the moment of inertia (/), the surface angular velocity (cj), and the surface rotational velocity (v rot ) of 
the companion star, respectively. 



ever, only a small fraction of the total angular momentum is dis- 
tributed over an extended outer layer (see Fig. IToV 

5. Comparison with Tycho G 

Tycho's SN (SN 1572) is one of only three historical SNe la 
observed in our galaxy. SN 1572 has the advantage that the field 
of SNR is not so crowded with stars and, therefore, it provides 
a good opportunity to observationally identify the companion 
in case of a SD progenitor. Here, we examine the viability of 
the candidate Tycho G as the possible surviving companion star 
in SN 1572 by comparing its observed rotat ional velocity of ~ 
6 ± 1.5 km s" 1 (IKerzendorf et alJl2009L [2012b with the results of 
our hydrodynamical simulations. 

In our MS_1 10 model, after the impact, the rotational veloc- 
ity of the star is significantly reduced to 23% of tha t bef ore the 
explosion of ~ llOkms" 1 (see Table[T]). In Section I47TI we dis- 
cussed that the surface spin of the companion star is roughly con- 
stant from 5000 s after the explosion. However, we stopped our 
simulation 2 x 1 4 s after the explosion because of the high com- 
putational costs. iPodsiadlowskil (120031) followed the post-impact 
evolution of the surviving companion star further and showed 
that the star might reesta blish thermal equ ilibrium ~ 10 3 yr af- 
ter the explosion (see also lPan et al.ll2012ah . Figure [9] shows that 
the radial distribution of the rotational velocity of the surviving 
companion becomes approximately constant ~ 1 x 10 4 s after 
the SN explosion, and the surface velocity at the equator is con- 



verged at the end of our hydrodynamics simulations. Therefore, 
we can safely assume that the surface velocity at the equator of 
the post-impact companion star would keep a stable value until 
its thermal equilibrium is established after a few thousand years. 
Considering that the Tycho SN remnant is only 439 years old, the 
spin of the MS_110 model (-25 kms" 1 ) after the impact is far 
larger than the observational ro tational velocity of the star Ty cho 
G (6 ± 1.5km s" 1 according to IKerzendorf et al]l2009i 12012) . 

iHanl (120081) carried out detailed binary evolution calcula- 
tions for the WD+MS channel of SNe la, in which RLOF starts 
when the companion star is on the MS or in the Hertzsprung-gap 
phase. They obtained many properties of the companion star at 
the moment of the SN explosion (e.g. their masses, spatial veloc- 
ity, effective temperature, luminosity, surface gravity, etc). These 
properties might be verified by the observations. The distribution 
of propertie s of companion stars in the plane of (v^, M^ N ) from 
IHanl (120081) is shown in Fig. I12h . Based on the observations of 
IKerzendorf et al.l (I2009ll2012h . the location of the rotation of the 
star Tycho G is shown with black error bars. Clearly, because of 
the bound rotation before the impact, Tycho G is far away from 
the allowed region. 

Next we estimate the rotational velocity of the stars shown 
in Fig. [T2li after the impact. In order to obtain the distribution 
of rotational velocities of surviving companions after the SN la 
explosion, we re-perform ID consistent binary evolution calcu- 
lations to construct the structures of the companions, obtaining 
additional three consistent models named "MS_081", "MS_131" 
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Table 1. SPH impact simulations for four different MS companion models. 



Model name 


Mf 
[M ] 


pSN 

[days] 


Rf 

[Ro] 


[Ro] 


V SN 
•rot 


<t 
[kms- 1 


vL 


rSN 

spin 

[10 50 | 


spin 

I cm 2 s _1 ] 


Abound 

[M ] 


MS_160 


1.21 


0.29 


0.93 


2.55 


160 


52 


98 


2.94 


1.31 


1.04 


MS_131 


1.23 


0.56 


1.45 


3.94 


131 


40 


78 


2.07 


0.92 


1.06 


MS_110 


1.18 


0.91 


1.97 


5.39 


110 


25 


46 


2.25 


0.62 


0.95 


MS_081 


1.09 


2.00 


3.19 


8.92 


81 


12 


16 


2.32 


0.26 


0.84 



Here, M 2 SN , P SN , Rf , a SN , vg? 



and / s s p ^ n are the mass, the orbital period, the radius, the spin velocity and angular momentum 



of the companion star at the moment of the explosion, respectively. vJ ot , / spin and M boU nd denote the spin velocity, the angular 
momentum, and the total bound mass of the companion star after the SN impact. vf ot is the rotational velocity at the surface 
after the thermal equilibrium is reestablished. Note that the rotational velocity, v^, is calculated by assuming that the rotation 
of the star is locked with the orbital motion due to tidal interactions. 
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Fig. 8. Left column: temporal evolution of the density distribution of the surviving companion star in model MS_1 10. Right column: 
temporal evolution of the radial configurations of the angular velocity, co (solid line), and the rotational velocity, v rot (dashed line). 
R is the distance from the rotation axes, the z-axis. Please note that we use different length scales in different diagrams. 
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Fig. 9. Radial profiles of angular velocity (solid curves) and 
rotational velocity (dashed curves) of the companion star 
in the MS_110 model at different times after the impact 
(0.0,0.5 •• -2.0 x 10 4 s). 



and "MS_160" with different rotational velocities of 81 kms" 1 , 
131 km s" 1 and 160 km s" 1 (see Table[T]). We then use these three 
models as input into our impact simulations to investigate the 
dependence of the rotational velocity of companion star after the 
impact on its value before the explosion. The properties of all 
companion models are shown in Table. [T] We find that to a good 
approximation the post-impact rotational velocity of star scales 
linearly with the pre-explosion velocity as is shown in Fig. [13] 
and can be fitted by 



V r f ot = 0.52 • V™ - 29.8 (kms" 1 ), 



(2) 



where V r f ot is the post-impact rotational velocity at the end of the 
simulations, and V™ corresponds to the star's initial rotational 
velocity at the moment of th e explosio n. 

In the previous work of iLiu et all (120121) . it was found that 
the unbound mass of the companion star caused by the SN im- 
pact is strongly dependent on the ratio of separation to the radius 
of the companion star, a^Ri. This relation can be fitted by a 
power law if the effect of the different structures of the compan- 
ion stars is neglected (s ee also lMarietta et al.l2000HPakmor et al.l 
120081: iPanet atll2012bl) : 



M , 



(3) 



where C\ is a fitting constant which depends on the different 
companion star models. The parameter J3 is the corresponding 
power-law index. 

Therefore, adopting above power-law relation (see also 
equation (2) of ILiu et al.l l2012h and the linear relation © ob- 
tained from the data shown in Fig. [13] we calculate the final 
bound mass and the post-impact rotational velocity of surviv- 
ing companion stars based o n all compan ion models from the 
population synthesis study of iHanl (120081) . The results are pre- 
sented in Fig. [T2b . The location of Tycho G is also shown with 
a black err or bars. Here, we assume th at Tycho G is a one solar 
mass star (iRuiz-Lapuente et all 120041) and its rotational veloc- 
ity is 6 ± 15 km s" 1 dKerzendorf e7al1l2012h . Fig. [T2l shows that 
Tycho G is located in the outer region of 95.4% of all systems, 
which casts doubt on Tycho G star as a promising candidate in 
SD progenitors of SNe la. 



Fig. 10. Radial rotation profiles of the MS_110 model 2 x 10 4 s 
after the explosion. R stal is the equatorial surface radius and v e is 
the equatorial rotational velocity. M T , J r , I r and oj y are the total 
bound mass, the total angular momentum, the angular velocity 
and the moment of inertia of all bound particles interior to the 
corresponding cylindrical surface, r, respectively. v r is the circu- 
lar velocity at distance r from the rotation axis. co c is the angular 
velocity on the rotation axes (z-axis). 




3.0 4.2 
r ( 10 11 cm ) 

Fig. 11. Evolution of the distribution of angular momentum in- 
side the companion star of model MS_1 10. The color scale gives 
the time since the explosion. Here, r is the distance from the 
rotation axes (z-axis), J v denotes the angular momentum of all 
bound particles interior to the corresponding cylindrical surface, 
r. Each colored solid curve shows the radial distribution of an- 
gular momentum, / r , at a given time. The dashed-dotted curve 
corresponds to the temporal evolution of the total angular mo- 
mentum of the companion star. 



However, it is not possible to exclude the star Tycho G com- 
pletely for several reasons. 

1. The errors shown in Fig. [12] are base d on an assumption 
that T ycho G is a one solar mass star (iRuiz-Lapuente et aTl 
120041) . However, it is very difficult to determine the actual 
mass from the observations. Therefore, Tycho G should be 
located in the vertical gray strip at its given rotation velocity 
in Fig. [T2b . If Tycho GisaO.6-0.7 M star rather a solar 
mass star, it could more likely be a candidate for the compan- 
ion star of SN 1572 in the SD scenario. Moreover, the errors 
given in Fig [12] are based on an assumption that the observed 
rotational velocity of Tycho G is reduced by an inclination 
angle of 60°. However, if the inclination angle is 30° instead 
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Fig. 12. Panel (a): distribution of rotation velocity versus mass 
of the co mpan i on sta rs at the moment of the explosion as ob- 
tained by iHanl (120081) . The number of systems decreases from 
the inner regions to outer ones. The blue area contains 50.0% of 
all stars, blue plus green 68.3%, and so on. Panel (b): as panel 
(a), but now the distribution is shown after the impact. The dis- 
tribution is shifted to lower rotational velocity and slightly lower 
masses. The unbound mass stri pped by the SN impact is cal- 
culated by using equation (2) of iLiu et all (120121) . and the post- 
impact rotational velocity of companion star is computed em- 
ploying the linear relation show n in Fig. [T3l( see t ext). The lo- 
cation of Tycho G according to iKerzendorf et al.1 (120091 12012b 
is shown with error in black assum ing a solar mass star with 
5% error (iRuiz-Lapuente et a"Dl2004 ) and an inclination angle of 
i = 60°. The black square gives the location of Tycho G for an 
assumed inclination angle of i = 30°. The red, vertical dashed 
line corresponds to the lowest rotational velocity covered by our 
hydrodynamical simulations (i.e., MS_081 model in Table[T]). 



(see black square in Fig. [T2b), the true rotational velocity 
of Tycho G would be higher, thus making the interpretation 
that it is the surviving companion of SN 1572 more likely. 
Figure [14] shows the distribution of the observable vj 0t sin i 
assuming a random viewing angle. 
2. In our simulations, we found that the post-impact rotation ve- 
locity of the surviving companions is linearly dependent on 
the initial rotational velocity at the moment of the explosion 
(see Fig. [13]). However, we restrict our simulations to cover 
a range of initial rotational velocities of only 80kms _1 < 
vf ^ < 160 km s" 1 (which corresponds to a range of the post- 
impact rotational velocities of lOkms" 1 < v f YOt < 50kms _1 ) 
because of computational costs. The companion star models 
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Fig. 13. Initial rotational velocity at the moment of the SN ex- 
plosion, V r s ^, versus rotational velocity after the impact, V r f ot , for 
four different MS companion star models (see Table [TJ. 

with lower initial rotation velocity always evolve to the gi- 
ant phase with larger orbital separation at the time when the 
supernova explodes, which significantly increases the com- 
putational effort of the simulations. FigfT2h shows that some 
companion- star models (~ 14% of all binaries) are located in 
a range of < 80 km s" 1 (i.e., regions left of the red verti- 
cal dashed line in Fig.[T2b). Actually, these models may not 
obey the linear relation as shown in Fig. [13] but rather fol- 
low something close to a power law relation. Therefore, the 
models with post-impact rotational velocities close to zero 
in Fig. [T2b might not be realistic. However, ~ 86% of all 
models are located in the range that we follow in the simula- 
tions (80 km s" 1 to 160 km s" 1 ), indicating that most of post- 
impact surviving companion stars in the WD+MS scenario 
should be rapidly rotating after the impact. 

3. The post-impact masses of the companion stars shown in 
Fig. [T2b are calculated by directly adopting the power law 
relation from equation (2) of ILiu et all (120121) . which ignores 
the effect of the different companion structures. However, 
ILiu et al.l (120121) found that the details of the companion 
structure also plays an important role in determining the 
amount of unbound mass caused by the impact. Nonetheless, 
the ratio of separation to the radius of companion star, <Zf /R2, 
is the most importan t parameter to determine the unbound 
mass (ILiu et al]l2012h . 

4. We assume that the surviving companion star has a con- 
stant surface rotational velocity after ~ 5000 s after the SN 
la impact. However, we only carried out our simulations to 
2 x 10 4 s after the explosion. Extended calculations based on 
the results of our impact simulations are still required to fol- 
low in detail the post-impact rotation of the star during its 
re-equi libration phase. For instance, very recently [Pan et al.1 
(l2012al) showed that the remnant star after the SN impact 
could continue to expand on a local thermal timescale of 
~ 10 3 yrs before it shrinks again. 

6. Rotation velocity after re-equilibration 

After the supernova impact, the companion star puffs up and the 
stellar envelope is out of thermal equilibrium. The post-impact 
remnant star will recover the state of thermal e quilibrium on the 
Kelvin-Helmholtz timescale of - 10 3 - 10 5 yrs (iPan et al.l2012at 
iMarietta etai1l2000t lPodsiadlow"skil 12003b . In order to estimate 
the rotation rate after thermal equilibrium is reestablished, we 
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Fig. 14. Same as Fig. [T2b , but showing the observable quantity 
v f rot sin i, where i is the inclination angle. The distribution of i 
is chosen corresponding to a random viewing angle. The black 
cross is again the location of Tycho G with an observed rotation 
velocity of 6 ± 1.5 km s" 1 . 



make the simplifying assumption that the remnant star has con- 
stant rotational velocity during the re-equilibration phase. After 
thermal equilibrium is established, the post-impact remnant star 
will shrink and return to the state of a MS -like star. As a result, 
the surface rotational velocity of the remnant star will increase 
again, and we assume that the angular momentum of the star 
would redistribute towards rigid-body rotation. If the angular 
momentum is given roughly by / = aMR 2 a> before the SN ex- 
plosion and after re-equilibration, and also has a constant value 
of the parameter a, the rotational velocity (V^ t ) of the remnant 
star after the re-equilibration can be determined once the radius 
of the re-equilibration star, R, is fixed. In order to estimate the 
final radius o f the remnant star, t he MS relation R oc M 2/3 is 
used (see also lMarietta et al.l200Qb . By adopting the post-impact 
angular momentum, J^ pin , and the value of a that is calculated 
for the pre-SN companion star model, we can estimate the rota- 
tional velocity of the remnant star after the thermal equilibrium 
is reestablished. The rotational velocity of different companion- 
star models before the explosion and after the re-equilibration 
can also be fitted with a linear relation in good approximation 
(see Fig. [T5b). In Fig.fTSb we also show the distribution of the 
rotational velocities, Vf oV after re-equilibration, employing the 
same method as in Section [5] It can be seen that the compan- 
ion stars relax to higher rotation velocity again after the thermal 
equilibrium is reestablished (at least about a few thousand years 
after the SN explosion). 



7. Summary and conclusions 

In this work, we have modeled the impact of SN la ejecta on 
their MS companion stars by means of 3D SPH hydrodynami- 
cal simulations which included the orbital and spin velocity of 
the companion stars. The MS companion stars were constructed 
by using Eggleton's ID stellar evolution code makin g use of the 
optically thick wind model of lHachisu et al.l (0~999). We found 
that the orbital motion and the spin of the companion star do 
not significantly affect the amount of unbound mass and the kick 
velocity caused by the SN impact. This result, obtained with a 
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Fig. 15. Top panel: similar to Fig. [13j but in the plane of (V^, 
Vf 0t ). Note that V r ^ t is the rotation velocity of the remnant star 
after thermal equilibrium is reestablished. Bottom panel: similar 
to Fig. Eb, but in the plane of ( V*, M\). 



SPH co de, differs from what was found previously by iPan et al.l 
who from their grid-based FLASH code simulations 
concluded that these two properties increase the amount of un- 
bound mass by up to ~ 16%, but also found that the kick veloc- 
ity is not affected. We have shown that the SN impact affects the 
rotation rate of the companion stars and their rotation laws. In 
our simulations we found that the SN impact removes as much 
as 55% to 89%, of the initial angular momentum due to the fact 
that 14% to 23% of the initial mass is stripped from the MS com- 
panion star. The remnant expands significantly after the impact 
which causes the equatorial surface rotational velocity to drop 
significantly to 14% to 32% of the original value. Additionally, 
we found that the post-impact rotational velocities of companion 
stars depend linearly on those prior to the explosion. Compared 
with the observed rotational velocity of the presumed companion 
star of Tycho's supernova, Tycho G, the final rotational velocity 
we obtain in our simulations is still higher by at least a factor of 
two. Whether or no this difference is significant, and may cast 
doubts on the suggestion that Tycho G is the companion of SN 
1572, has to be investigated in future studies. In particular, hav- 
ing a more accurate mass and inclination of Tycho G would help, 
as well as simulations that follow the evolution for much longer 
than we can do with our explicit SPH c ode. Further more, by us- 
ing the population synthesis results of iHanl (120081) . we showed 
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the distributions of the rotational velocities of surviving com- 
panions after the impact, which may be very useful for further 
observations to identify the surviving companion stars in SNRs. 
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